Negative Mobility induced by Colored Thermal Fluctuations 
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Anomalous transport of non-Markovian, thermal Brownian particle dynamics in spatially-periodic 
symmetric systems that is driven by time-periodic symmetric driving and constant bias is investi- 
gated numerically. The Brownian dynamics is modeled by a Generalized Langevin equation with 
exponentially correlated Gaussian thermal noise, obeying the fluctuation-dissipation theorem. We 
study the role of non-zero correlation time of thermal fluctuations for the occurrence of absolute 
negative (linear) mobility (ANM) near zero bias, negative-valued, nonlinear mobility (NNM) and 
negative differential mobility (NDM) at finite bias away from equilibrium. We detect that a non- 
zero thermal correlation time can either enhance or also diminish the value of ANM. Moreover, 
finite thermal noise correlation can induce NDM and NNM in regions of parameter space for which 
such ANM- and NNM-behavior is distinctly absent for limiting white thermal noise. In parts of the 
parameter space, we find a complex structure of regions of linear and nonlinear negative mobility: 
islands and tongues which emerge and vanish under parameters manipulation. While certain such 
anomalous transport regimes fade away with increasing temperature some specific regions inter- 
estingly remain rather robust. Outside those regimes with anomalous mobility, the ac/dc driven 
transport is either normal or the driven Brownian particles are not transported at all. 

PACS numbers: 05.60.-k, 05.45.-a, 85.25.Cp, 74.25.Fy 



I. INTRODUCTION 



Out-of-equilibrium transport processes are of promi- 
nent interest in modern statistical physics from all, fun- 
damental, experimental and device aspects pi, !2]. A typ- 
ical far from equilibrium problem is generally nonlinear 
response of a system to an external stimulus: For exam- 
ple when an external constant force Fq acts on the system 
and all other forces are of zero-average, it is expected that 
the long-time, stationary average particle velocity v be- 
comes an increasing function of the load Fq. An everyday 
example is that of a force pushing objects on a table. A 
linear Ohmic resistor characteristic constitutes another 
example: An increase of voltage is accompanied by a lin- 
ear increase of current. This normal response behavior 
is distinct from cases with anomalous transport. Famil- 
iar examples are the emergence of negative differential 
mobility or conductivity, or the nonlinear response in- 
volving negative (-valued) nonlinear mobility away from 
the linear response behavior around zero applied volt- 
age. Here our focus is on yet another anomalous trans- 
port behavior, namely so called "absolute negative mo- 
bility" . This latter anomalous transport behavior refers 
to a regime where the resulting velocity or current as- 
sumes the opposite sign of the applied force or voltage 
around the zero bias regime. While the NDM is common 
for an abundance of nonlinear systems, the phenomenon 
of ANM has been experimentally detected predominantly 
much less frequently. Some examples that come to mind 
are the nonlinear response in p-modulation-doped GaAs 
quantum wells , or also semiconductor superlattices [4J , 
and recently also in driven Josephson junctions [5]. The 
phenomenon of ANM can typically relate to a genuine 



quantum effect, involving asymmetry tunneling dynam- 
ics. In contrast, ANM as a result of classical stochastic 
dynamics occurs more rare; but it is expected to occur 
whenever stylized, ratchet-like structures including geo- 
metric entropic barriers are present. A flurry of recent 
theoretical [HI [7] and less few experimental works E] 
indeed prove such behavior. 

The effect of ANM can occur also in the form of a far 
from equilibrium phenomenon in driven nonlinear sys- 
tems such as in nonlinear underdamped Brownian motion 
dynamics [9j I10|, 111] or even in overdamped nonlinear 
Brownian motion in presence of time-delayed feedback 

m- 

In this work we shall focus on the case of time- 
dependent driven underdamped Brownian motion occur- 
ring in a periodic, reflection symmetric potential and 
driven by thermal correlated noise. In prior works [51 110] . 
we have studied the transport properties of a classical 
Brownian particle of mass m moving in a spatially peri- 
odic potential V{x) = V{x + L) of period L and barrier 
height AV , which is subjected to an external unbiased 
time-periodic force F{t) = F{t + T) of period T = 27r/f] 
with angular frequency f2 and of amplitude A. Addi- 
tionally, a constant bias Fq acts on the system. This so 
deflned Brownian particle dynamics is then modeled by 
the driven Langevin equation [S]; i.e., 

mx -f 7i = — V'{x) 

+ Acos(m + 0o) + ^o+ CW, (1) 

where x = x{t) is a position of the particle at time t, 
a dot denotes differentiation with respect to time and 
a prime denotes a differentiation with respect to the 
Brownian particle coordinate x. The parameter 7 de- 
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notes the viscous friction strength and (f>Q is an initial 
phase of the time-periodic driving. Here, the thermal 
fluctuations are modeled by zero-mean, Gaussian white 
noise ^(i) with the Dirac delta auto-correlation function 
{£,{t)£,{s)) = 2kBToS{t — s), where ks the Boltzmann con- 
stant and To denotes the temperature. 

We could show that in the above system, there are 
distinct regimes of anomalous transport. In particular 
we can identify both (i) ANM and (ii) NDM and the 
phenomenon of (iii) negative nonlinear mobility (NNM), 
occurring away from the linear response regime with re- 
spect to external bias Fq. We also remind the reader 
that ^ mimics the physical realization for the behav- 
ior of a physical Josephson junction |9j [TOj E] . In this 
latter case, the periodic potential V{x) has the explicit 
sinusoidal form, i.e., 

V{x) = AVsm{2nx/L). (2) 

Notably, the theoretical findings for ANM in Josephson 
junctions has recently been verified experimentally with 
the work in Ref. [5]. 

II. PERIODICALLY DRIVEN AND BIASED 
NON-MARKOVIAN BROWNIAN DYNAMICS 

The thermal noise ^(t) in Eq. ([ij is approximated to 
be ideally white noise with zero noise correlation time 
Tc — 0. In real systems, however, the correlation time 
of thermal fluctuations is only approximately zero. This 
approximation is justified if Tc is much smaller than the 
smallest characteristic time of the system itself. There 
are many examples where this situation is well satis- 
fied in real systems. There are, however, also situa- 
tions where the thermal correlation time Tc is of order or 
greater than r^, so that the white-noise approximation 
fails [131 [131 [TS]. In this latter case a modeling based 
on the Markovian Langevin equation ([T|) is not correct; 
instead, the generalized Langevin equation should then 
be invoked. 



A. The Generalized Langevin dynamics 

When the thermal noise is correlated, the appropriate 
Langevin dynamics is a non-Markovian dynamics with 
memory-friction described by the so called Generalized 
Langevin Equation (GLE) [TOJ [T71 [M [HI [20] . It explic- 
itly reads 

mx{t)+ [ K{t-s)x{s)ds^-U'{x{t),t)+^{t), (3) 
Jo 

where the full potential takes the form 

U{x, t) ^V{x)~ [A cos(m + </>o) + Fo]x. (4) 

This non-Markovian dynamics can be derived from first 
principles by means of coupling the system of interest 



to a bath of harmonic oscillators [TTJ [TO] with the total 
system being prepared in canonical thermal equilibrium 
It then follows from the central limit theorem that 
the thermal fiuctuations ^(i) obey a zero-mean, station- 
ary typically non-Markovian Gaussian stochastic process. 
The auto-correlation function of the thermal noise f (t) is 
related to the memory (frictional) kernel K(t) via the 
fluctuation-dissipation relation 16, 17, 1S[ [TOllMj : 

mm) = kBToK{\t^s\). (5) 

Interestingly, due to the nonlinearity of the potential 
V{x) it is still an unsolved, open problem to derive the 
explicit form of the generalized master equation for the 
single-event non-Markovian probability p(a;, see in 
Refs. [TSj, [HI [2T|; this task is achieved only in form 
of a time-convolutionless master equation with time- 
dependent transport coeflicients iff the potential V{x) 
is at most quadratic in x only jl9l I21j : then yielding a 
general non-Markovian Gauss process for the equilibrium 
dynamics {x(t), x(t)). 

B. Exponentially correlated thermal fluctuations 

The Gaussian thermal fluctuations ^(i) in ^ are com- 
pletely determined by the memory function K(t). If 
the memory function K(t) is the Dirac delta function, 
i.e., K{t) = 2jS{t), then Eq. ^ reduces to the form 
(jlj. A well studied form of correlated fluctuations is 
defined by means of an Ornstein-Uhlenbeck (0-U) sta- 
tionary stochastic process for ^(t) [2]. This Gaussian 
Markov process for the thermal noise (note that the re- 
sulting Brownian dynamics is then still non-Markovian) 
is henceforth correlated exponentially. We next use a 
Markovian embedding of the GLE dynamics in Eq. (|3| . 
Towards this objective we present the correlation func- 
tion in the form 

mm) = ksTom - s\) = e-l*-^l/^% (6) 

Tc 

where is the correlation time of the 0-U process which 
we can then smoothly vary from the limit of white Gaus- 
sian noise (tc — 0) to strongly correlated thermal noise 
(tc ^ To). Because the integral kernel exhibits an ex- 
ponential form, we can convert Eq. (|3| into a set of 
ordinary stochastic differential equations: Let us define 
the auxiliary stochastic process w{t) via the relation 

w{t) = ^ / e-(*-^)/"=i(s) ds. (7) 

Tc Jo 

Then Eq. ([3| is equivalently transformed into the form 

mv{t) = ~U'{x{t),t) ~w{t) + C{t), (8) 
x{t) = v{t), (9) 

w{t) = --w{t) + ^v{t), (10) 

Tc Tc 

m = --m + -V^^ik^oT{t), (11) 

To Tc 
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where the normaHzed white Gaussian noise T{t) obeys 
(r(t)r(s)) = S{t — s) while the last equation of this set 
describes the 0-U noise with the exponential correlation 
function (|6| [21 [13]. Note that in Eq. the linear 
combination z{t) = £,{t) — occurs. By subtracting 
the two last relation we then find the set of three cou- 
pled Markovian Langevin equations, modeling the two- 
dimensional non-Markovian GLE in Eq. (Isl; i.e., 



x{t) 
v{t) 

m 



v{t), 

1 



U'{x{t),t) + -z{t), 



(12) 
(13) 



-z{t) - ^v{t) + T{t). (14) 



The corresponding three-dimensional Fokkcr-Planck 
equation for p{x, v, z, t) is numerically cumbersome to im- 
plement for V{x) a nonlinear, spatially periodic function. 
Alternatively we apply direct numerical methods for the 
solution of the three coupled Langevin equations in Eqs 



The limiting case of white noise has been analyzed in 
detail in Refs [HI , where a numerical method has been 
described, see also in Refs ^24j i25j. In the remaining of 
this work we shall use the dimensionless form of the set 
of equations ( 12 1-( 14 1. In doing so we scale coordinate x 



and time t as follows 



t 

t ^ — . 

To ' 



Then, the set of equations (12)-(14l is recast as 



X 
Y 

Z 



= Y 



-W'{X) +acos{Lui+ 0o) + / + Z, 



-Z 



-Y 



where 



Y 



L 



AV 



(15) 

(16) 
(17) 

(18) 
(19) 



and a dot denotes differentiation with respect to the re- 
scaled time i. Here, Y — Y{i) is the dimensionless veloc- 
ity of the Brownian particle and Z — Z{t) denotes the 
corresponding dimensionless random force. The remain- 
ing re-scaled parameters are: (1) the friction coefficient 
7 = (7/™)'''o = tq/tl equals the ratio of two character- 
istic times, namely time tq and the relaxation time of 
the velocity degree of freedom, i.e. tl = m/j; (2) the 
potential W{X) = V{x)/AV = W{X + 1) = sin(27rX) 
possesses unit period and barrier height AV — 2; (3) the 
amplitude a — LA/AV and the frequency to = ^tq (or 
the period T = 27r/w); (4) the load / = LFq/AV; (5) the 
zero-mean white noise ^{t) is correlated as {£,{t)£^{s)) = 
6{t — s) with a re-scaled noise intensity D — kBTg/AV. 
The latter is given as the ratio of two energies, the ther- 
mal energy and the half of barrier height of the potential 
V{x). 



From here on, we shall use only these dimensionless 
variables and shall omit the notation "hat" in all quanti- 



ties in Eqs (16l-(18l 
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FIG. 1; (color online) Regime of destructive role of finite ther- 
mal noise correlations on absolute negative mobility (ANM): 
The averaged long-time velocity v of the ac driven Brownian 
particle is depicted as a function of the externally applied, 
static force /. The chosen system parameters are: a = 4.578, 
uj = 4.9, 7 = 0.9 and D = 0.001. While white thermal noise 
(tc = 0) does depict ANM this is increasingly diminished 
upon increasing the correlation time of thermal noise until it 
is undone into a normal, positive mobility behavior. 



III. NUMERICAL RESULTS FOR BROWNIAN 
NON-MARKOVIAN ANOMALOUS TRANSPORT 

We next study numerically the long-time transport 
characteristics of the nonlinear Brownian dynamics with 
memory friction. In particular we shall focus on the cur- 
rent feature as given by the long-time averaged velocity 
V = (Y). This averaging is performed as follows: First 
we perform an average over all realizations of the ther- 
mal fluctuations which yields a temporally varying quan- 
tity. A second temporal average is the cycle period of the 
external ac-driving. Because the resulting asymptotic, 
long time dynamics is not necessarily ergodic (i.e. inde- 
pendent of chosen initial conditions) in all phase space 
we also need to perform an average over unbiased initial 
conditions. We have chosen uniformly distributed initial 
positions X{t = 0) = Xq over one period of the peri- 
odic potential W{X); the initial velocities Y{t = 0) — Yq 
are unbiased and taken as uniformly distributed in the 
interval [—2,2] and the initial phase 0o G [0j27r]. The 
symmetry consideration of Eq. ([3| then implies that this 
average velocity v(f) as a function of the external con- 
stant force / is an odd function, i.e., v{—f) = —v{f); 
thus v(f = 0) = 0. Therefore we will in our numerics 
only consider the half- axis with / > 0. 

Because, as shown with previous works with white 
thermal noise [Hllinilll]) the transport dynamics becomes 
very rich indeed in all parameter space exhibiting all, 
namely ANM, NDM and NNM. As it must be expected 
this richness does not diminish with yet another param- 
eter of variation, namely the correlation time of thermal 
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noise Tc- Also it must be kept in mind that it is im- 
possible to scan numerically over all possible parameter 
space of the driven nonlinear Brownian dynamics. In the 
following we shall focus our numerical study to regimes 
in parameter space that are (i) experimentally accessible 
[9l Ho] and (ii) exhibit a most interesting complexity for 
the nonlinear response. 

Clearly, the resulting velocity v — v{f) is typically non- 
linear in external bias /. The linear response behavior is 
defined for small bias / — > as 



vif ^ 0) = ^lf, 



(20) 



where the (linear) mobility fi can become negative, /i < 0. 
This regime will be termed absolute negative mobility 
(ANM); i.e. an anomalous transport regime for which the 
the particle is transported in the opposite direction to the 
externally applied force /. Moreover, negative nonlinear 
mobility (NNM) refers to an anomalous transport regime 
for which we find w(/)// < in some finite intervals of /, 
being disjoint from the interval around / = . Finally, 
regimes of negative differential mobility (NDM), when 
dv{f)/df < in some intervals of /, can be detected. 



A. Controlling ANM with noise correlation time 

1. Undoing ANM with small noise correlation time 

We first consider small bias values / and study the 
linear (/ — > 0) response and the accompanying nonlinear 
response with increasing / to larger values for specific pa- 
rameter settings. We first zoom into a parameter regime 
for which we find ANM for Tc = 0. An example is de- 
picted with Fig. 1. If the correlation time increases, 
starting out from zero, we observe a diminishing ANM 
with increasing Tc, until it disappears and turns into nor- 
mal, positive- valued mobility upon increasing Tq further. 
Put differently, the mobility coefficient /i in Eq. (20) 



starts to increase from negative values, passes through 
zero and eventually becomes positive. For Tc > 0.035, 
the velocity monotonically increases in the region of small 
bias / with > 0. Note that the dimensionless corre- 
lation time is rather short in comparison to the other 
characteristic time scales, that is with the characteris- 
tic time-scale tq = 1, the characteristic velocity relax- 
ation time-scale is I.Itq and time of periodic driving is 
r = 1.28ro. 

A first main finding therefore is that even a small, fi- 
nite thermal noise correlation time can diminish and even 
undo ANM. 



2. Enhancing ANM and creating NNM 

Upon scanning the parameter set to a different forcing 
strength of the ac-driving we find that thermal noise cor- 
relations can in fact also enhance rather than diminish 
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FIG. 2: (color online) Regime of constructive role of thermal 
noise correlations on absolute negative mobility: The aver- 
aged long-time velocity v of the Brownian particle is depicted 
as a function of the bias force / for selected values of the 
noise correlation time Tc of thermal fluctuations. The sys- 
tem parameters are chosen different from Fig. [4] and read: 
a = 4.1293, uj = 4.9, 7 = 0.9 and D = 0.001. Panel (a): 
Within this regime we find that a small thermal noise cor- 
relation time can enhance ANM. Panel (b): Upon further 
increasing thermal noise correlation time Tc from the values 
shown in (a) we observe a turnover into ANM again, note 
the case with Tc = 0.07; followed by a NNM behavior, see 
Tc = 0.0985, and finally a re-entrance into a normal transport 
regime is observed at small bias with a regime following with 
NDM as / is increased further; note the case with Tc = 0.12. 



the value of ANM. This is illustrated in Fig. |2] panel 
(a), for noise correlation times varying in the interval 
Tc e [0,0.028]. Beyond = 0.028 ANM diminishes 
again, exhibiting as well regimes with NNM, cf. the case 
= 0.05. 

It is intriguing to note that within this parameter set- 
ting a further increase of the correlation time yields an 
opposite behavior, see panel (b) in Fig. 2. The mobility 
coefficient \i then is decreasing from positive values and 
next turns into ANM again. The most pronounced ANM- 
value occurs around around Tc = 0.08 (not shown) before 
entering a regime with coexistence with both ANM and 
NNM around Tc — 0.1. For larger correlation time we 
find normal linear mobility followed up in the nonlinear 
regime with a region exhibiting NDM. 

We note that the various situations described above do 
occur and can coexist in other settings of the parameters. 
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FIG. 3: (color online) Regions of negative (nonlinear) mobil- 
ity manifested either as ANM or also NNM, in the param- 
eter plane given by {a, Tc} for three different temperatures, 
D = 0.0001,0.001,0.01 (top to bottom). The remaining pa- 
rameters are fixed: namely a = 4.1293, to = 4.9, 7 = 0.9 and 
/ = 0.1. 



B. Correlation time induced Islands and Tongues 
of ANM and NNM 



While in the previous subsection we have varied the 
external bias / away from zero we next keep this bias 
fixed at / = 0.1 and scan instead the correlation time Tc 
versus the ac-driving strength a. The emerging asymp- 
totic, averaged velocities are then depicted in color-coded 
plots as shown with Fig. [3] for three different settings of 
the temperature Tq. 



Not unexpected, there occur in this parameter space 
several domains with anomalous transport features. The 
geometric structure of these domains in the depicted 
regime of the {a, Tc}- variation is very complex. Let us 
remind ourselves that the underlying deterministic dy- 
namics is chaotic in some regimes and therefore a fractal 
structures of certain domains must be expected to ex- 
ist. We are interested in the stability of those domains 
in parameter space for which ANM and/or NNM occur. 

Our comprehensive numerical investigations reveal a 
rich diversity of structures, formed by regions of both, 
ANM and NNM. We present one example and visualize 
it in the parameter-plane {a, Tc}. In Fig. 3, we depict 
how the plane is divided into regions of normal {v{f) > 
for / > ) and anomalous (w(/) < for / > 0) transport. 
In doing so we here do not discriminate between ANM 
and NNM. Both these transport behaviors are jointly pre- 
sented. Two thermal situations, namely, 'low' and 'high' 
temperatures are shown for comparison. At 'low' tem- 
perature, cf. panel (a), we reveal the refined structure 
with many narrow, slim and twisted regions of ANM- 
NNM. Some of those regions, the 'tongues', survive with 
correlation time Tc approaching zero (i.e. the horizon- 
tal abscissa) and there are 'islands' of negative mobility 
which disappear for = 0. If one fixes one of the pa- 
rameters, say Tc = 0.1, intervals of negative mobility are 
clearly noticeable: there are several intervals of the am- 
phtude a G (0^,0^+1) for which ANM-NNM can be de- 
tected. Outside these intervals, a normal response to the 
load / is found. 

If temperature is increased, this small temperature 
structure is increasingly washed out; i.e. it becomes 
more smooth. Many previously existing domains of 
ANM-NNM behavior start to shrink or vanish altogether. 
We detect some few robust regimes for which anoma- 
lous transport persists, namely a few islands and a few 
tongues. The most robust such island against increasing 
noise strength (temperature) is the island located around 
{a, Tc} — {6.56, 0.31}. This very domain even survives at 
'high' temperatures. 

We emphasize that such complicated regimes of ANM- 
NNM are not just rare occurrences: they can be veri- 
fied with numerically arbitrarily-high-accuracy calcula- 
tions and over extended intervals in parameter space. 
Given the complexity of the underlying dynamics with 
time-dependent ac-driving, nonlinearity and in presence 
of noise with finite correlation time the observed behavior 
is clearly beyond a sensible analytical description. 



C. Current-voltage characteristics of a realistic 
Josephson junction device 

As an application of the above theoretical study we 
consider next a Josephson junction for which the anoma- 
lous conductance have been measured in Ref. 0. The 
relation in Eq. ([l| with the potential ([2| models the resis- 
tively and capacitively shunted Josephson junction, also 
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known as the Stewart-McCumber model [1511211155]. It 
contains three additive current contributions: a Cooper 
pair tunnel current characterized by the critical current 
/q, a normal (Ohmic) current characterized by the nor- 
mal state resistance R and a displacement current due to 
the capacitance C of the junction. For this model, the 
position X of the Brownian particle translates into is the 
phase difference (j) between the macroscopic wave func- 
tions of the Cooper pairs on both sides of the junction, 
i.e., X — (f>, the mass m = {h/2e)'^C, the friction coeffi- 
cient 7 = {h/2e)^{l/R), the barrier hight AV = {h/2e)lQ 
and the period L — 27t. The load Fq — {h/2e)Id is 
given by means of the dc-bias current, the amplitude 
A — {h/2e)Ia and the frequency define the external 
ac current. The velocity v = x translates into the volt- 
age across the junction. The experimental results for 
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FIG. 4: (color online) Role of non-zero thermal noise cor- 
relations on absolute negative mobility for a parameter set 
corresponding to the anomalous Josephson junction transport 
regime used in the experiment in Ref. [5J . The long-time aver- 
aged voltage V across the junction is depicted as as a function 
of the externally applied dc current /. The corresponding di- 
mensionless parameters are: a = 13.874, ui = 2.75, 7 = 2.264 
and D = 0.000895. The white thermal noise regime (i.e. 
Tc = 0) depicts ANM-behavior which is sustained for short 
correlation times of the thermal noise. For increasingly larger 
correlation times ANM is diminished and eventually vanishes 
all together and crosses over into a normal, positive mobility 
behavior. Still, negative-valued, nonlinear mobility can oc- 
cur far away from equilibrium; note the nonlinear response 
behavior for r = 0.02 and around / ~ 0.6. 

the Josephson junction presented in Ref. |5] show very 
good agreement with the Stewart-McCumber model with 
the following set of parameter values: namely, the di- 
mensionless amplitude of the ac-current a = 13.874, the 
frequency ui = 2.75, the friction coefficient 7 — 2.264 
and the (white) noise intensity D — 0.000895. We 
note that Eq. ^ can serve as a generalization of the 
Stewart-McCumber model towards a regime of validity 
to lower temperatures where finite correlations of the 
thermal fluctuations increasingly play a signiflcant role 
|28j . We numerically find that at small non-zero correla- 
tion time of the thermal fiuctuations ANM is sustained 



within tailored parameter regimes, as studied previously 
in the context of current-voltage characteristics of a ac 
driven Josephson junction device in Ref. |S], see Fig. 
4. With increasing thermal noise correlation time ANM 
is weakened and finally turns into normal behavior fol- 
lowed by a regime of NNM far away from equilibrium. 
This feature of finite thermal noise correlation is of rel- 
evance from an experimental point of view when mea- 
suring mobility as a function of temperature. We can 
also conclude from this study that at the temperature of 
the experiment performed in Ref. [5], the white thermal 
noise approximation is seemingly well satisfied. 



IV. CONCLUSIONS 

With this study we numerically analyzed the role 
of non-zero correlation time of thermal fluctuations on 
the anomalous transport regimes of underdamped, non- 
Markovian Brownian particles that are driven by time- 
periodic and static forces. We detected a rich variety of 
anomalous transport behavior in an experimentally wide 
parameter space where anomalous transport can be mon- 
itored. The regions of absolute negative mobility and 
negative nonlinear mobility form complicated structures 
in parameter space with stripes, fibres and islands, see 
Fig. 3. At low temperatures these structures can unam- 
biguously be attributed to finite correlations of thermal 
fluctuations. Such correlations are either constructive or 
destructive in nature with respect to the size of anoma- 
lous transport coefficients such as mobility. A subsequent 
increase of temperature starts to blur these structure. 
Nevertheless, there occur stable 'islands' with anoma- 
lous, negative-valued mobility behavior which are solely 
induced by non-zero thermal noise correlations, cf. Fig. 
3. 

In the regime of linear mobility response, such non- 
zero thermal correlations are found either to diminish or 
also to enhance the value of absolute negative (Hnear) 
mobility, dependent on the specific parameter setting for 
ac-driving strength and /or remaining parameters. 

We also have compared our predictions in a parame- 
ter regime of a recent experiment on anomalous response 
behavior in a ac/dc driven Josephson junction [5^: The 
observed ANM behavior is sustained for small thermal 
noise correlation time but increasingly fades out with in- 
creasing thermal noise correlation, see in Fig. 4. 
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